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By using the recent nonperturbative equation of state of the quark gluon plasma 
{SJ . derived in the formalism of the Field Correlator Method, we investigate the bulk 

' properties of the strange quark matter in beta-equilibrium and with charge neutrality 

at T = p = 0. The results show that the stability of strange quark matter with 
■ respect to 56 F e is strongly dependent on the model parameters, namely, the gluon 
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condensate G2 and the qq interaction potential V%. A remarkable result is that the 



^ width of the stability window decreases as V± increases, being maximum at V± = 

q; and nearly zero at V\ = 0.5 GeV. For V\ in the range < V\ < 0.5 GeV, all values 

CD 4 

of G2 are lower than 0.006 — 0.007 GeV obtained from comparison with lattice 



results at T c (/u = 0) ~ 170 MeV. These results do not favor the possibilities for the 
existence of (either nonnegative or negative) absolutely stable strange quark matter. 
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I. INTRODUCTION 

Since the works of A. R. Bodmer [l| and E. Witten j^j, the existence of strange quark 
matter (SQM) has been largely investigated. SQM is a particular form of matter comprising 
roughly equal amounts of u, d and s quarks. It is presumed that SQM has been produced 
at extreme conditions of high temperatures and densities in the beginning of the universe 
and/or latter at low temperatures and high densities in compact stellar interiors (e. g., 
neutron stars and quark stars). In this regard, it is of importance to mention the pioneer 



work o: 



a. 



N. Itoh, who investigated the possibilities for the existence of hypothetical quark 



stars 

According to the Bodmer- Witten conjecture, SQM might be more stable than ordinary 
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nuclear matter. In a pioneer work, E. Farhi and R. L. Jaffe J4] investigated SQM in equilib- 
rium with respect to weak interactions, at zero temperature and pressure, in the context of 
the MIT Bag Model p]. In this model, quarks enter in the respective equation of state (EOS) 
as a free quark gas with a Fermi-Dirac distribution. The confinement is represented by a 
bag enclosing the free quarks with a constant B which gives the vacuum energy density dif- 
ference between the confined and deconfined phases. Improvement of the model is obtained 
by the inclusion of corrections to first order in the QCD coupling constant in the a c < 1 
(perturbative) regime (see [4| and references therein). Until now, SQM properties in dense 



nuclear matter and compact stars interiors have mostly been considered in the 



xamework 



of the MIT Bag Model \4. 16l-l20|. Applications of the Nambu-Jona-Lasinio 2l|, [22[ quantum 



field theoretical approach have also been done to describe quark matter properties in com- 



pact stars interiors 
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-|27|. In alternative investigations, in which quark masses are 



density dependent, SQM is also taken as a free Fermi gas mixture of quarks and anti-quarks 
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291 ] . These approaches are used to describe u, d and s quark matter at zero and (not 
so high) nonzero temperatures and large density regions where the approximation of free 
quarks can be considered. However, this is not so at all densities and (low) temperatures. 
Quarks strongly interact subjected to a potential that is large when the qq distances are 
large. In this case, nonperturbative methods must be considered. 

QCD, the fundamental theory of strong interactions, due to its nonlinearity has been 
taken as an inappropriate theory for the purposes of practical applications as, for instance, 
the calculation of an EOS to describe quark matter at all finite temperatures and densities, 
including nonperturbative effects of confinement. Asymptotically, the description of quark 
matter becomes simple, but at low temperatures and (or moderate) densities the attempts 
to obtain an EOS including the confinement have appeared as a difficult task to be achieved. 
However, since 1987, the investigations of Yu. A. Simonov [30] and H. G. Dosch 31 33 1 
have resulted in the construction of an important method based on vacuum field correlators 
functions that is being continuously developed up to now. 

Recently the nonperturbative EOS of quark-gluon plasma was derived in the framework 
of the Field Correlator Method (FCM) |3J|, also called Stochastic Vacuum Model. FCM 
(for a review see 35| and references therein) is a nonperturbative approach which naturally 
includes from first principles the dynamics of confinement in terms of color electric and color 



magnetic correlators. The parameters of the model are the gluon condensate G2 and the 



3 



qq interaction potential V\ which govern the behavior of the EOS, at fixed quark masses 
and temperature. FCM has been used to describe the quark-gluon plasma dynamics and 
phase transition 36l-l39t|. An important feature of the model is that it covers the entire phase 
diagram plane, from the large temperature and small density region to the small temperature 
and large density region. In the connection between FCM and lattice simulations, the critical 
temperature at // c = turns out to be T c ~ 170 MeV for G 2 ~ 0.006 - 0.007 GeV 4 Q, Q. 
In a previous work, the existence of stable SQM on strange star surfaces (see Sec. 3. 3 in 

n 

[40l |) was shortly considered in the FCM framework. In the present work, we perform a 
more detailed investigation of SQM on the same line. The system is a gas of u, d and s 
quarks and gluons subjected to the interaction potential V\. The vacuum energy density 
difference between confined and deconfined phases, A|e mc |, is a nonperturbative quantity 
expressed in terms of G 2 . The main parameters of our calculation are A|e mc |(or G 2 ) and Vi, 
and the strange quark mass m s (assuming m u = = 0). We first investigate the energy 
per baryon, from which we obtain the stability window with respect to the 56 F e nucleus. 
Strangeness, hadronic electric charge and density are also addressed. 

This application of the FCM to the study of the bulk properties of SQM (not considered 
before) shows the role of the method to provide alternative indications for its parameters. 
In 40j we have shown the importance of the comparison of the FCM calculations with 



astrophysical observations of some strange star candidates. Similarly, in the present work, 
the main purpose is the relevance of the FCM predictions for the SQM properties to be 
compared with lattice simulations and/or the results at RHIC and LHC experiments. 

This paper is organized as follows. In Sec. [Ill we summarize the theoretical framework 
of the FCM and show the equations to be used in our calculation. In Sec. Ill II we show the 
results and in Sec. IIVI we give the final remarks and conclusions. 



II. BASICS EQUATIONS 



re confined-deconfined phase transition is dominated by the 



In the FCM approach, t 
nonperturbative correlators 35(. The dynamic of deconfinement is described by Gaussian 
(quadratic in F^ U F^ U ) colorelectric and colormagnetic gauge invariant Fields Correlators 
D E (x), Df(x), D (x), and Df{x). The main quantity which governs the nonperturbative 
dynamics of deconfinement is given by the two point functions (after a decomposition is 
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made) 



g 2 { ir f [Ei(x)$(x, y)E k (x)<$>{y, x)] 
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where u = x — y and 



y) = Pexp 



(2) 



(3) 



is the parallel transporter (Schwinger line) to assure gauge invariance. 

In the confined phase (below T c ), D E (x) is responsible for confinement with string ten- 
sion a E = (1/2) J D E (x)d 2 x. Above T c (deconfined phase), D E (x) vanishes while Df (x) 
remains nonzero being responsible (toghether with the magnetic part due to D H (x) and 

(x)) for nonperturbative dynamics of the deconfined phase. In lattice calculations, the 



nonperturbative part of Df(x) is parametrized in the form 35l. |41| 



Df(x) = Df (0)e 



-MM 



(4) 



where A = 0.34 fm (full QCD) is the correlation length, with the normalization fixed at 
T = v = 0, 



D E (0) + Df(0) 



where G2 is the gluon condensate [421 ]. 
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(5) 



The generalization of the FCM at finite T and \x provides expressions for the thermo- 
dynamics quantities where the leading contribution is given by the interaction of the single 
quark and gluon lines with the vacuum (called single line approximation (SLA)). As in [401 ]. 



from 



34( and standard thermodynamical relations 43j, we here explicitly rewrite in more 



convenient forms (for our purposes) the expressions (for one quark system, Nf = 1) for the 
pressure 



P, 



SLA _ 1 2A C 
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energy density 
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and include the number density of the quark system 



n 
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L SLA (T,» q ,J?) 



djJL q (27r) J 
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d A k 



f q SLA (T,J?,v q ) + f q SLA (T,jF, N ) 
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e P(E+TJ?-» q ) + i 



and f q LA (T, Hq., jf 



1 



e /3(E+TJf + x 



(9) 



(g = u,d,s), E = ^Jk 2 + m|, /3 = 1/T, and Jf = V\/2T is the exponent of the Polyakov 
loop in the fundamental representation, where 



pp poo 

Vi= / dr(l-rT) / td£D?(y/? 
Jo Jo 



(10) 



is the large distance static qq potential 
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4l|. 



The pressure and energy density of gluons are given by 



and 



-SLA 



SLA 
Pgl 



3Pgl 



V. 



(27T) 



d A k- 



k 



e /3(fc+Tjf ) _ i 



T 2 ^(iV c 2 -l)- 



d 3 k 



e /3(fc+TJf ) _ i 



(11) 



(12) 



fj-f is the exponent of the Polyakov loop in the adjoint representation. In 



91 * <9T v ^' c ^(2tt) 3 

where = | 

Eqs.©, ( fill and ( fl2l) . when Vi = we recover the ordinary Fermi and Bose gases. In order 
to give Eqs.(l6T)- (Tr2"|) in its most general forms, it was assumed that V\ is, in principle, a 
function of temperature and chemical potential. However, according to the parametrization 
given by Eq. ffl"0]) . V\ does not depend on the chemical potential. As pointed out in jsjj, the 
expected /i-dependence of V\ should be weak for values of // much smaller than the scale of 
vacuum fields (which is of the order of ~ 1.5 GeV) and is partially supported by the lattice 
simulations 44j. As in 40|, |45|, |46], we take V\ independent of the chemical potential, so 
dJ?/dn q = in Eqs. (GO) and ©. 

In order to take into account the presence of electrons to keep the quark matter in /3- 
equilibrium and with charge neutrality, we also include the equations for the pressure, energy 
density and number density of electrons given by 



1 2 
3(2^ 



d 3 k— [/ e (T,/x e ) + / e (T,/x e )] 
hj 



(13) 
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where 



e e = J d 3 k E e [f c (T,fx e ) + / e (T,/ie)] , (14) 

He = (2ffp / ^M/eC^e) - 7e(T,/i e )] , (15) 



/ e (r, // e ) = e/3(ge _ Me) + l , / e (T, / X e ) = e/3(gc+ | tc) + i (16) 



and E c = + . 

The composition of SQM is maintained in /3-equilibrium with respect to weak interactions 
and in electric charge neutrality. The weak interactions reactions are given by 

d — > u + e + v e (17) 

and 

s^u + e + v e . (18) 

As pointed out in 4|, the neutrino gas is so dilute that it play no role in the dynamics 
of the system. So, by neglecting the neutrino chemical potential, the chemical equilibrium 
equations are given by 

/id = /i« + /ie (19) 

and 

fi s = [id ■ (20) 
The overall charge neutrality requires that 

l(2n s u LA -n s d LA -n s s LA )-n e = 0. (21) 

By numerically solving Eqs. ( IT9|) -( l2T|) . for each value of the input total density 

n = n SLA + n SLA + n SLA + n ^ (22) 

the unknown chemical potentials fi u , /id, yU s and /i e are determined for fixed values of T, G 2 



401 ]. as explained 



and V\. However, our calculation here is slightly different from that in 
below. 

The total pressure and energy density of the quark-gluon system, including electrons are 
given by 

P = Pf A + E P"^ ' ^ac\ + Pe , (23) 
q=u,d,s 
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.SLA 



.SLA 



where 



Ale, 



q=u,d,s 
11 



3*> 



-AG 2 , 
32 



(24) 



(25) 



is the vacuum energy density defference between confined and deconfined phases in terms 
of the respective difference between the values of the gluon condensate, AG 2 = G 2 (T < 



T c ) — G2(T > T c ) ~ \Gi [36|, [37|, and Nf is the number of flavors. 

We follow the same line of 4j to investigate the behavior of the SQM at zero temperature 1 
and total pressure by solving the above equations for constant values of the energy per 
baryon, E/A = e/n^, where = (n u + + n s )/3 is the baryon number density. We 
perform our calculation in the m u = = approximation 2 (in Sec. IHH larger values of m u 
and rrid are speculated in order to look for strangeness excess and negative electric charge 
possibilities). By using the constant E/A constraint in the above equations, m s and A|e„ ac | 
(or G 2 ) are determined for fixed values of V\ in the range < V\ < 0.5 GeV (where V± = 0.5 



GeV is the value of V% at T = T c obtained from lattice investigations |47|). As in [40 
our calculation here is made for V\ constant, so T 2 dJf /dT = —V±/2 in Eq. ([7j) 



45 



46|, 



A. Quark matter at T = and constant Vi. 



For pedagogical purposes, we show the previous equations for the quark system at zero 
temperature and constant V\ for the general case of nonzero quark masses. Zero temperature 
implies that 

f^iT.^J?) _^ e(ti q -E-TJ?) 

T-»0 

and Eqs.©-® lead to 



(26) 



P. 



SLA 



N c k 3 q / , 2 , 2 3 2 



k q Jk% + m* -m 2 g In 



kg + \Jk 2 + 



(27) 



1 In reality, we take T = 0.001 GeV in Eqs. (|6l TT5l) which is a good approximation or, alternatively, by using 
Eqs.([2l| - 

2 For our purposes here, it is irrelevant if electrons are assumed massless or not. 
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k q + ^Jkq + 



2\ 1 



m„ 



n 



SLA 



3 
JZ 
2 3 ' 



iV, fc. 

7T : 



K = \] - v i/ 2 ) 2 - m \ , (9 = u,d,s). 
When Vi = 0, the ordinary Fermi momentum fc^ is recovered. 



(28) 



(29) 



(30) 



1. Zero mass approximation 



In order to better understand the role of G2 and V% in the study of stability of quark 
matter, it is instructive to apply the above equations to quark matter at zero temperature 
and pressure, assuming that all quark species are massless particles. This simple case serve 
to help us to understand the behavior of the constant E/A curves in Fig. (TJ as well as the 
shrinking of the stability window in panel (a) of Fig. [31 at m s = 0. For massless quarks, 
Eqs.(E7D-([29]) (with iV c = 3), by using Eq.fl3D}, are reduced to 

1 



4tt 2 



Vx 



(31) 



.SLA 



1 f3 



Vi 



and 



7T 2 1 9 2 



(32) 



(33) 



At total zero pressure, the sum of the quark pressures is balanced by the vacuum energy 
density, 



1 ( Vl 



(34) 
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and the energy density is 

e = y^ j e q + A\e v 



3 E"« + |5E(^-t) S+ ^ 



q 
3 

= - Vt n A + AA\e vac \ (35) 

where (here, in this section) the baryon number density is n A = {n u + n<i + n s )/3 for SQM 
and tia = {n u + n^)/3 for nonstrange quark matter. Notice the presence of the additional 
term (3/2) Vin^ with respect to the corresponding expressions in the MIT Bag Model 
Also noticed is that the sum of the quark pressures and energy densities are given in terms 
of A|£ mc | (or G2) and V\ (differently from the MIT Bag Model where they are given solely 
in terms of the bag constant B). Now, let us particularize the above equations for two and 
three flavor quark matter. 

Two flavor 

For a gas of u and d quarks, charge neutrality (neglecting the not important cotribution 
of electrons as in [9j) requires that n s d LA = 2 n^ LA , from which it follows that (/^ — V\/2) = 
2 1 ^ 3 (Hu — Vi/2) . The two-flavor vacuum energy density (for Nf = 2 in Eq. (l25l) ) is A|e TOC | U( i = 
(29/192) G 2 - Thus, the energy per baryon of the ud system is 

= (1 + 2 4 / 3 ) 3 / 4 (4vr 2 ) 1/4 ( A\e vac \ ud )*/< + \ Vl 

A J ud 1 

= 6Ml{A\e vac \ vd ) 1/ * + ~V 1 

= 4.016 G 2 1/4 + -Vi . (36) 



Three flavor 

The three flavor quark system (SQM) is naturally charge neutral, with n^ LA = n^ LA 



n 



SLA 



t*u = (J>d — Us = A* > an d n e — A*e = . The vacuum energy density (for Nf = 3) is 
A|e mc | = (9/64) G 2 , and the energy per baryon becomes (using the previous notation for 
the SQM system, without subscripts) 



(§) = s^Y^a^i)^^ 

= 5.714 ( A\e vac \ + Ivi 



2 

3.499 G 2 1/4 + ^Vt . (37) 
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In Eqs. (l36l) and (157)1 . for fixed E/A, the increase of V(G 2 ) is compensated by the corre- 
sponding decrease of G 2 (Vi). As shown below in Fig. [TJ (at m s = and a given E/A) the 
maximum value of G2 is obtained for V\ — . 

The energy per baryon of 56 F e is 930.4 MeV, so in this simple analysis the stability 
of SQM relative to iron corresponds to G 2 < (0.266 — 0.428V) 4 . As a result, we obtain 
G 2 < 0.005 GeV 4 for V l = 0, G 2 < 0.002 GeV 4 for V x = 0.1 GeV, and so on, until a very 
small value of G 2 10~ 5 GeV 4 ) for V = 0.5 GeV. This behavior explains the shrinking of 
the stability window at m s = shown in panel (a) of Figj3j 

Finally, from Eqs. ( 136]) and ( 1371) we obtain 

(E/A) 3 3/4 (47r2) l/4 ( A | £me | } l/4 + L5 ^ 

(E/A) ud (1 + 2 4 /3)3/4( 47r 2)i/4( A\e vac \ ud )V4 + 1.5 Vl 



5.714 ( A 




| ) 4 / 4 + 1.5 Vi 


6.441 ( A| 




ud )V*+ 1.5 V 



= 3 ' 499 G +L5l/ ' . (38) 
4.016 G 2 1/4 + 1.5 Vi 

It is evident that (E/A) < (E/A) ud (for the same values of G 2 and V\ in SQM and ud 
systems). From the last line of Eq.flHE]) it follows that (E / A) / (E / A) ud = 0.87 for V = 0. On 
the other hand, for the MIT Bag Model, with the correspondences A|e,; ac | = A\e vac \ ud = B 
and V = 0, we obtain (E / A) / (E / A) ud = 0.89 as in Q. 



III. RESULTS 



We are concerned with the bulk properties of SQM and concentrate ourselves to investi- 
gate the stability with respect to the 56 F e nucleus. In our investigation, m s enters as input 
parameter and A|£ mc |(or G 2 ) is determined for fixed values of E/A and V\. By this way, as 



in 



40j , but with a different logic, we obtain a scenario for the model parameters, indepen- 



dently of what the results of lattice calculation may be. We discuss the relations between 



the parameter values required 
with lattice predictions in jsj], 



or the SQM stability and the values obtained by comparison 



37]. 



In Fig. CEJ the constant E/A contours give m s vs A|e„ ac | (for the purpose of comparison 
with MIT Bag Model results in {4, 8, 9|) for different values of V. In order to understand the 
role the gluon condensate, we use the relation between A|£„ ac | and G 2 , given by Eq. (125]) . to 
plot m s vs G 2 for the same values of V. The contours are very sensitive to the values of V, 



11 



being shifted towards lower values of A|e vac | and/or G 2 as shown for V\ = (panels (a) and 
(b)) and V 1 = 0.01 GeV (panels (c) and (d)). To the right of the E/A = 0.93 GeV contour 
(in reality, 930.4 MeV corresponding to the energy per nucleon of 56 F e ), SQM is unstable 
with respect to the iron nuclei. The vertical line at the left of each panel is the limit of 
A|e V ac| and/or G2 when m s becomes large, so the strangeness per baryon goes to zero (see 
panel (a) in Fig. [2]). In this case, there is no distinction between strange and non-strange 
quark matter. Contours with E/A < 930.4 MeV terminate at the crossing with the vertical 
line of 56 F e . For V\ = 0, the results shown in panel (a) are numerically equivalent to the 
ones found in the case of the MIT Bag Model. However, we remark that in the FCM the 
vacuum energy difference A|£ vac | is essentially a nonperturbative quantity given in terms of 
the gluon condensate. Also shown is the E/A = 0.939 GeV contour corresponding to the 
nucleon mass. 

Stability window of the SQM is the region of allowed values of m s and A|e vac | (or G 2 ) 
where the energy per particle is lower than the one of 56 F e (bounded by the E/A = 0.93 
GeV contour and the respective vertical line). In the FCM, the stability of SQM depends 
on the values of V\ and/or G2. For a given value of E/A, the higher V\, the lower G2 (cf. 
Eq. (I37p for the case m s = 0). Moreover, even for V\ — (for which the contours present 
the maximum G2 at m s = 0), the values of G2 within the stability window are lower than 
0.006 — 0.007 GeV 4 obtained from lattice data on the critical temperature (sg), l^]. The 
possibility of the SQM be more bound than 56 F e is realized only for G2 < 0.005 GeV 4 . This 
has also been the case for m u = 5 MeV, = 7 MeV and m s = 150 MeV for which we have 
shown that G2 < 0.0041 GeV 4 for the existence of stable SQM in strange star surfaces 40]. 
Even if we take stability with respect to the nucleon mass (E/A = 0.939 GeV), instead of 
56 .F e , the values of G2 remain lower than the one in 0, \&\ . 



Fig. [2] shows for the given value of E/A (the same as in |4|) the strangeness per baryon 
(defined as in Q]) in panel (a) as function of m s for some values V\. The strangeness per 
baryon is always lower then unity, going to zero at m s ~ 0.3 GeV (however, strangeness 
excess might be possible for nonzero m u and/or as shown in panel (a) of Fig. [5]). For 
the same values of Vi, panel (b) shows the decrease of the baryon number density from 
its maximum at m s = until a constant value around m s ~ 0.3 GeV. We must be aware 
that for larger values of V\ and at some value of m s , the baryon number density might 
becomes lower than a critical value (if it exists) at which the phase transition takes place. 
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However, the determination of such a critical value is not the scope of the present work. 
We have also calculated the hadronic electric charge per baryon, 



shown in panel (c). For m s = the equilibrium configuration is given by an equal number 
of u, d and s quarks (n u = rid = n s and n c = 0) with zero electric charge. When m s and/or 
V\ grow, the system develops a positive hadronic electric charge. For large m s , the hadronic 
electric charge per baryon saturates at a constant value which also depends on V% (cf. panel 
(b) of Fig. H]). For V\ = 0, this saturation point is ~ 0.0056 at m s ~ 0.3 GeV as in 4]. 

Given that all E/A < 0.93 GeV contours are within the stability window, it is instructive 
to consider some features of SQM at E/A = 0.93 GeV of 56 F e . In panel (a) of Fig. [3] we 
show that the overall effect of the confining forces is to shift the stability windows towards 
lower values of G2 for increasing values of V\. The windows not only narrow as m s grows 
at a fixed V\, but they also narrow as V\ grows at the same value of m s . In particular, for 
V\ = 0.5 GeV the vertical line is located at a negligible value of G2 and the stability window 
width is very small (not visible in the scale of the figure) . From the locations (at fixed values 
of Vi) of the constant E/A contours and the respective vertical lines on the horizontal axis 
at m s = 0, where each window presents its maximum width, we construct two plots V\ vs 
G2 as shown in panel (b). The region between the dashed and solid curves illustrates the 
decrease of the stability window width with V\. In panel (c), we show the baryon number 
density for < V\ < 0.5 GeV. For V\ = 0.5 GeV, it is nearly zero. So, as we have observed 
above, due to the decreasing of it might happens that a quark-hadron phase transition 
occurs at some value of V\ and m s . 

Fig. H] shows the strangeness per baryon (in panel (a)) and the hadronic electric charge 
per baryon (in panel (b)), at the energy per baryon of 56 F e , as function of m s for V\ — 
and V\ = 0.5 GeV. For all values of V\ in the region < V\ < 0.5 GeV and m s in the region 
< itl s ^ 0.35 GeV, the strangeness per baryon is always less than unity. Depending on the 
values of V\, the saturation of the hadronic electric charge per baryon is between 0.0056 for 
V\ — and ~ 0.006 for V\ = 0.5 GeV (this variation is not visible in panel (c) of Fig. [2] 
because of the low values of V\). 

Let us now consider the question of the strangeness excess. We have performed our 
calculation assuming that m u = = 0. However, for nonzero m u and/or m^, strangeness 



Z 2n u -n d -n 



s 




(39) 



A n u + n d + n 
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excess can be obtained which is more sensitive to rrid than it is to m u . For the usual 
values of u and d quark masses, m u = 5 MeV and rrid = 7 MeV, the strangeness excess 
is less than 1%, but it can be larger for larger m u and m^. As an ilustrative example, we 
have (speculatively) extrapolated m u and rrid beyond its usual values in order to obtain a 
strangeness excess around 9-14 % as shown in panel (a) of Fig. EJ for the given quark masses, 
V\ and E/A. Generally speaking, this excess only occurs for low values of m s (<, 0.02 GeV) 
and large values of V\ which, in turn, correspond to very small values of G 2 (~ 10" 6 GeV 4 ) 
(the increase of u and d quark masses also shifts the E/A contours towards lower values of 
G 2 ). 

Correspondingly, we also speculate the possibility of negative hadronic electric charge. 
The change of the hadronic electric charge is more sensitive to m u and m s than it is to 
rrid- For large values of V±, it can happen that strange quarks are more abundant than the 
massless u and d quarks in the region of small values of m s . So, negative hadronic electric 
charge appears to be allowed for large values of V\, m u and rrid (as for the case of strangeness 
excess), but for small m s , as shown in panel (b). In this case, instead of electrons, a sea 
of positrons neutralizes the negative hadronic electric charge. We also remark that, at the 
same values of Vi, m u and ma, larger values of S/A and lower (negative) Z/A are allowed 
for E/A < 0.899 GeV and m s in the region < m s < 0.02 GeV. 

Summarizing the above results, strangeness excess and negative hadronic electric charge 
per baryon are realized only for large values of rrw and/or rrid and V\ (say, Vitt 0.3 GeV), 
but for values of G 2 much lower than the one in jsfl, Ei^ . For m u = rrid = as well as for the 
usual m u = 5 MeV and rrid = 7 MeV, we observe that the values of the model parameters 
obtained in the present paper do not favor the existence of neither nonnegative nor negative 
SQM with energy per baryon lower than the one of the 56 F e nucleus. 



IV. FINAL REMARKS AND CONCLUSIONS 



In this work we have investigated the bulk properties of SQM by using the quark-gluon 



plasma EOS derived in the FCM nonperturbative approach [34J. The important parameters 
of the model are the gluon condensate G 2 (which enters the EOS through the vacuum 
energy difference A|e vac | between confined and deconfined phases) and the large distance qq 
interaction potential V\. The results have shown that confinement plays an important role 
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for the stability of SQM. 

We have performed the calculation in the m„ = = approximation and assumed 
SQM in /3-equilibrium and charge neutrality, at zero temperature and pressure. We have 
been mainly concerned with the absolute stability of SQM with respect to 56 F e nucleus. In 
order to look for stability windows of SQM, we have started our investigation by drawing 
contours of constant energy per baryon, E/A, in the m s vs A|e vac | (and/or m s vs G 2 ) plane 
for fixed values of V\. Strangeness and hadronic electric charge have also been considered. 
Our study revealed remarkable features which we summarize as follows. 

The general trend is that the SQM stability is very sensitive to the values of the model 
parameters responsible by the confining forces. A remarkable aspect is that the behavior 
of the stability window strongly depends on the values of V\. For increasing values of this 
parameter, the constant E/A contours and also the respective stability windows as a whole 
are shifted towards lower and lower values of A|e: vac | (and/or G 2 ) until to ~ at V\ — 0.5 
GeV. Moreover, the width of the stability windows diminish when V\ becomes larger. At 
m s = 0, it has the maximum width between G 2 = 0.003 GeV 4 and G 2 = 0.005 GeV 4 
for Vi = and a nearly zero width at G 2 ~ for V\ = 0.5 GeV. This amounts to say 
that absolutely stable SQM would exists, in principle, for < V\ < 0.5 GeV (although 
somewhat problematic at V\ = 0.5 GeV due to the smallness of the corresponding value 
of G 2 ). However, a striking point is that the values of G 2 are lower than the one in the 
range 0.006 - 0.007 GeV 4 obtained from comparison with the lattice data at the critical 



temperature [36|, |37{|. This point puts a severe restriction for the existence of absolutely 
stable SQM with respect to 56 F e . 

We have also calculated the strangeness per baryon and the hadronic electric charge per 
baryon. As m s grows, the strangeness per baryon decreases from S/A = 1 at m s = to 
S/A = for some value of m s which depends on E/A and V\. Correspondingly, the hadronic 
electric charge per baryon is always nonnegative, rising from Z/A = up to a constant 
value which depends on the value of V\. Another remarkable feature (in the m u = rrid = 
approximation) is that S/A < 1 and Z/A > for all values of m s within the stability 
window. 

In the attempt to find strangeness excess and negative hadronic electric charge, we have 
observed that S/A > 1 and Z/A < appear to be allowed only for very large values of 
m u and/or (beyond the usual ones), but for small m s (<, 0.02 GeV) and large V\ (say, 
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between 0.3 GeV and 0.5 GeV). However, the corresponding values of Gi remain lower than 
the one in as in the case of m u = = 0. 

From the above, in the context of the FCM approach, it appears that the values of the 
model parameters obtained in our investigation do not favor the existence of absolutely 
stable SQM. Of course the above results depend on the constant V\ assumption, in the 
present work. 

Taking into account the importance of experiments at RHIC and LHC, it is instructive at 
this point to consider finite temperature effects on the SQM properties. In order to check the 
influence of nonzero temperatures on the SQM stability at zero pressure, we have applied 
the same procedure employed for the T = case to study the behavior of the stability 
window shown in Fig. [3j but for T ^ 0. We have taken several values of T up to 30 MeV 
for constant V\ and for V\(T) parametrized in [37] for T >T C as 

for T = T c , 2T C , 3T C and some arbitrary values of T c along the phase diagram transition 

n 

curve. In both cases, the results are qualitatively analogous to those in Fig. 1 of [8[ and Fig. 
3 of |29|. However, V\(T) is decreasing with the growth of T, so the shift of the stability 
window towards lower values of G2 takes place as T — > T c . For V\(T = T c ) the result is the 
same as for constant V\ = 0.5 GeV at T = shown in panel (a) of Fig. [3j 

Generally speaking, our results appear to be consistent with the fact that absolutely 
stable SQM has not been observed up to the present. The experiment with STAR at RHIC 
has not confirmed the existence of SQM nor proved that it does not exist 48-5(3]. The low 



values of G2 with respect to the ones in 36|, |37J] would provide a possible explanation for 
the absence of absolutely stable SQM signature. However, before any conclusion towards 
the nonexistence of absolutely stable SQM, we must have in mind that our theoretical 
results should be a consequence of the approximations contained in the development of the 
FCM nonperturbative EOS. FCM is a robust theoretical approach where the dynamics of 
confinement is one of the most important aspects of the model. Therefore, we must be 
aware that the FCM nonperturbative EOS is presently developed in the so called Single 
Line Approximation, where the confinement dynamics include only single quarks and gluons 
interactions with the vacuum (34]]. 

On the other hand, in our calculations, V\ and G 2 were taken as /i-independent pa- 
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rameters. As pointed out in |45|, |46J, the /i-independence of Vi should be a questionable 



assumption. Also, in the large density domain, important effects should be related to a 
possible density dependence of G2 [5l|- In our opinion, these aspects are very interesting 
possibilities to be considered. However, this does not have been the scope of the present 
work. 
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FIGURE CAPTION 

Fig. 1 - Contours of constant E/A of strange quark matter at zero temperature and 
pressure for two different choices of the parameter V\ {E/A and V\ in GeV units). The 
vertical dashed line is where the energy per baryon number of the two-flavor quark matter 
exceeds the one of 56 F e . Panels (a) and (c): strange quark mass as function of Ae vac . Panels 
(b) and (d): strange quark mass as function of The contours labeled (N) stand for the 
nucleon mass. The 0.93 contours correspond to 930.4 MeV energy per nucleon of 56 F e . 

Fig. 2 - Panel (a): The strangeness per baryon as function of the strange-quark mass 
for different choices of V\. Panel (b): As in panel (a), but for the baryon number density. 
Panel (c): As in panel (b), but for the hadronic electric charge per baryon. 

Fig. 3 - Panel (a): Stability windows of SQM at zero temperature and pressure, bounded 
by the E/A = 0.93 GeV contour of 56 F e and its vertical line, for different choices of V± (in 
GeV units) labeling each contour. Panel (b): for each value of V\ in panel (a), the respective 
locations of the E/A contour (solid) and its vertical line (dashed) on the horizontal axis. 
Panel (c): Baryon number densities at the energy per baryon of 56 F e as function of the 
strange quark mass for different values of V\. For V% = 0.09 GeV and V% = 0.167 GeV, we 
have Ua = no at m s ~ 0.3 GeV and m s = 0, respectively, where Uq = 0.153 fm -3 is the 
nuclear saturation number density. 

Fig. 4 - Panel (a): Strangeness per baryon as function of the strange-quark mass for the 
given values of V x (in GeV units), all at E/A = 0.9304 GeV of 56 F e . Panel (b): As in panel 
(a), but for the hadronic electric charge per baryon. 

Fig. 5 - Panel (a): Strangeness per baryon as function of the strange-quark mass for 
the given nonzero u and d quark masses, V\ and E/A (all in GeV units). The energy per 
baryon E/A = 0.899 GeV was taken for the purpose of comparison with the results in |4|. 
Panel (b): As in panel (a), but for the hadronic electric charge per baryon. 
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